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The effects of twin boundaries (TBs) on the complex interaction between magnetism and su- 
perconductivity in slightly electron-doped Ba(Ca)(FeAs)2 superconductors are investigated. The 
spatial distributions of the magnetic, superconducting and charge density orders near two different 
types of TBs are calculated. We find that TBs corresponding to a 90° lattice rotation in the a-b 
plane enable magnetic domain walls to form with only a small effective Coulomb interaction between 
valance electrons, and that superconductivity is enhanced at such TBs. Contrastingly, we find that 
superconductivity is suppressed at TBs corresponding to an asymmetrical placement of As atoms 
with respect to the Fe atoms in the a-b plane. 
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The recent discovery of Fe-pnictide based superconduc- 
tors offers an alternative avenue to explore the physics of 
high temperature superconductors [iHSl- Similar to the 
cuprates, the parent compounds of the Fe As-based su- 
perconductors also possess antiferromagntic (AF) ground 
states 0, With increasing electron or hole dop- 

ing, the AF order is suppressed and superconductivity 
(SC) appears in both the cuprates and the Fe-pnictides. 
However, different from the cuprates, SC and a 2 x 1 
collinear AF or spin- density- wave order can coexist in 
doped Ba(FeAs)2 superconductors 0, 0|- Because each 
unit cell of these new materials contains two inequivalent 
Fe ions, different organizations of the magnetic moments 
of Fe ions in both normal and superconducting states can 
lead to a diverse assortment of magnetic structures and 
unusual electronic properties jl, @| . 

Recently, twin boundaries (TBs) oriented 45° from 
the x(a)-axis were observed in the normal state of 
Ca(Fei_ x Co :c As)2 [10]. Across these TBs, the a-axis of 
the crystal rotates by 90°, and the modulation direction 
of AF order that exists is rotated by 90° as well. That is, 
90° magnetic domain walls (DWs) are formed at the TBs. 
Also, in the SC state of underdoped Ba(Fei_ ;E Co ;E As)2 
with x < 0.07, it has been found that the diamagnetic 
susceptibility is increased and that the superfluid density 
is enhanced on the same type of TB [ll|. C onsistent with 
these experiments, a theoretical study [12| found that 90° 
DWs can be formed at low doping levels and that SC is 
enhanced on them. However, the DWs considered in that 
study formed in the absence of TBs, being induced in- 
stead by a strong effective Coulomb interaction between 
valance electrons, while in the experiments fiol fill ] the 
DWs were pinned at TBs. In this letter, in order to better 
understand the effects that TBs have on the development 
of magnetic and superconducting order, we investigate 
the magnetic, SC and charge density orders near 4 dif- 
ferent TB configurations using the Bogoliubov-de-Gennes 
(BdG) equations for very under-doped Ca(or Ba)(FeAs)2 
compounds. Among these TBs are ones involving assy- 



metric placement of the As atoms above and below the 
Fe plane. This type of TB has not been previously stud- 
ied. This study is based upon the band model [13| and 
the phase diagram for electron doped Ba(Fei_ ;E Co ;E As)2 
[14j . We predict that the enhancement or suppression of 
SC, the location of DWs and the electron-density distri- 
butions are largely dependent on the nature of TBs. Our 
results provide a theoretical explanation for the relation- 
ship between TBs and enhanced SC order observed in 
experiments. 

Consider a Hamiltonian H = Hq + Hsc + Hi nt that 
describes the energy of valance electrons. Ho is a 
non-interacting energy from a two-orbital tight-banding 
model, the detailed form of which can be found in 
Ref. 13, T3 |. The pairing interaction energy of the elec- 
trons is 

H SC = ]T A^c^cj^ + H.c. 
i/ij no- 
where A^jj, is the pairing parameter between two elec- 
trons, one at position i with the orbital \i and the other 
at position j with orbital and cj^ is the creation oper- 
ator of an electron with spin a at position i with orbital 
li. Here a denotes the opposite spin of a. The mean-field 
magnetic interaction energy 
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H int = (U-3J H ) 



where U is the on-site Coulomb interaction, Jh is the 
Hund's coupling, n\^ G is the electron number operator, 
and (riipa) is the local electron density. The eigenvalues 
and eigenfunctions of the total Hamiltonian H can be 
obtained by self-consistently solving the BdG equations 
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where 

H- ll2 j U(7 = -tirfv + [U {nipa) + (U - 2J H )(m fl a) 

+(U - 3J H )(nij2a) - tolS-^S^ 

is the matrix element of H with the same spin a between 
the orbital \i at position i and the orbital v at position j, 
and to is the chemical potential. The pairing parameter 
Ai^jjy and the local electron densities (n^) and (n^) 
satisfy the following self-consistent conditions 
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where V\^ v is the pairing strength and f(x) 



SC order pa- 



is the 

Fermi-Dirac distribution function. The 
rameter at position i is defined as Ai 
Ai 5 i_£_£Ai 5 i + £_£+Ai 5 i_£ + £), the local magnetic moment 
at position i is defined as m\ = |X^(( n i/it) — ( n i/4))? 
and the total charge density at position i is defined as 
(m) = X^(( n i/it) + The chemical potential to is 

determined by the electron filling per site n (n = 2 + x), 
and for the value of the hopping terms t\^ v are assumed 
to be ti_ 4 = 1,0.4,-2,0.04 HEQ- Only the electron 
pairings of the same orbital between the next-nearest- 
neighbor Fe sites are considered. For example, we choose 
Vinjv = 1-4 for ii = v and |i — j| = V% and zero for all 
other cases. This choice of the pairing potential implies 
that the SC order has s±-wave symmetry 15, 16| . 

The phase diagram of the electron-doped 
Ba(Fei_ :c Co :c As)2 compounds as a function of tem- 
perature T and doping x has been qualitatively mapped 
out with U = 3.4 and J H = 1.3 Q. When both T and 
x are small, 2x1 collinear AF order is found, with SC 
order uniformly distributed over the sample. However, 
this AF order is unstable against the formation of the 
90° magnetic DWs oriented 45° from the x-axis as the 
strength of U is increased toU = 4.8 for small x at T = 
(l2| . In the following, the spatial profiles of magnetic, 
SC and charge density order near four different types 
of TB at T = will be investigated. Throughout this 
work, we set x = 0.04, U = 3.8 and Jh = 1.3. Note that 
DWs do not form spontaneously in the absence of TBs 
at these parameter values, instead the superconducting 
and charge density orders are uniform and AF order 
exists. 

A. Diagonal TB of the Lattice- The lattice in the 
Fe-plane of these compounds is almost square, having 
slightly different lattice constants a and b along x- and 
^-directions [TtJ • TBs can thus be formed by exchanging 
the lattice constants a and b on the opposite side of the 
TB. Figure 1(a) shows the structure of a single such TB 
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Figure 1: (a) The lattice structure near a diagonal TB (red 
dashed line), the open circles represent the positions of Fe 
atoms, a (blue solid line) and b (black solid line) are the 
lattice constants along x and y directions in domain A. Cor- 
responding spatial profiles of (b) the magnetic order, (c) the 
superconducting order, and (d) the charge density order are 
presented. 



oriented at 45° with respect to the x-axis. Since there is 
only a small difference in the magnitudes of a and 6, this 
TB can be accounted for by assuming slightly different 
nearest neighbor hoping terms t a = 1.0 and t^ = 1.2. 
To analyze the effects of this type of TB, we considered 
a 28 x 28 lattice with periodic boundary conditions di- 
vided into 4 different domains separated by three parallel 
TBs (not shown) along the lines y = x + 14, y = x and 
y = x — 14. As shown in Fig. 1(b), there are three 90° 
DWs formed and pinned on the TBs. The patterns of 



these quantities are very similar to those found in 12] 
without the TBs (see Fig. 2 (a) to 2(c) in [12]). However, 
in this case, the DW forms with a smaller value of the 
effective Coulomb interaction U, indicating that the ex- 
istence of this type of TB is beneficial to the formation 
of the 90° DWs. The solutions presented in Fig.l are al- 
ways stable against the uniform 2x1 collinear AF order 
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Similar to the results without TBs [12], the SC as well 
as the charge density get significantly enhanced on the 
DWs, which occur at the TBs, and suppressed in the 
middle of the magnetic domains (see Fig. 1(c) and 1(d) V 
All of this is in good agreement with experiments 10, Hi- 
lt is important to note that the lattices on both sides of 
the TB should be well matched at the TB, and that each 
of the unit cells along the TB is only slightly deformed 
from the square shape. Thus, we do not expect that 
scattering of the electrons from any disorder due to the 
TB would be strong. 
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Figure 2: (a) The lattice structure near a TB parallel to the 
^/-direction (red dashed line), the open circles represent the 
positions of Fe atoms, a (blue solid line) and b (black solid 
line) are the lattice constants along x and y directions in do- 
main A. Corresponding spatial profiles of (b) the magnetic 
order, (c) the superconducting order, and (d) the charge den- 
sity order are presented. 
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Figure 3: (a) The lattice structure near a twin boundary (red 
dashed line) parallel to the direction formed by misplac- 
ing As atoms. The open circles represent the positions of 
Fe atoms, and the red and blue dots respectively denote the 
As(up) and As(down) atoms. Corresponding spatial profiles 
of (b) the magnetic order, (c) the superconducting order, and 
(d) the charge density order are presented. 



B. Parallel TB of the Lattice- A TB formed by ex- 
changing the a and b lattice constants can also be ori- 
ented parallel to the x- or ?/-axis (as shown in Fig. 2 (a)). 
We studied this case by considering a a 28 x 28 lattice 
with periodic boundary conditions divided into 3 differ- 
ent domains separated by two TBs (not shown) along the 
lines x = 7 and x = 20. 

Here the magnetic DWs are pinned at the TBs (see 
Fig.2(b)) on which weak local ferromagnetic order ap- 
pears. The SC has a periodic modulation and is en- 
hanced on the DWs, but suppressed in the middle of the 
magnetic domains (see Fig.2(c)). A charge density wave 
appears near the DWs (see Fig.2(d)) while the electron 
density gets suppressed in the middle of the magnetic 
domains. It is important to point out that in this case, 
the lattices on the opposite sides of a TB are not well 
matched. Therefore, there may be considerable scatter- 
ing of the electrons due to the disorder near these TBs. 
If this effect is included, we expect that the SC would get 
suppressed, instead of being enhanced, on the DWs and 
the TBs [ill. 

C. Parallel twin boundary due to the asymmetry of As 
atoms- Another possible TB can be generated by slip- 
ping the lattice on the right side of the TB by a lattice 
constant along the y-direction with respect on the lat- 
tice on the left of the TB. There are two different types 
of As atoms in our model, we label them as As (up) and 
As (down) atoms relative to the Fe plane. This can be 
clearly seen from Fig. 3(a), in which the TB is represented 



by the red-dashed line. The crystal lattice for the FeAs 
layer has T>2d symmetry, namely the 4 nearest neighbor- 
ing As atoms of a "down" As atom should be all "up". 
The hopping terms between the next-nearest-neighboring 
Fe ions via the hybridization of the 4p orbital with the 
As atom in the middle should have different values de- 
pending on whether the As atom is above (£2) or below 
(£3) the Fe plane [l3|, EH- The D2d symmetry is broken 
by the presence of the TB. We considered a 28 x 28 lat- 
tice with periodic boundary conditions divided into three 
domains by two TBs located at x = 7 and x = 20 (not 
shown). The lattice constants are both a along x- and 
y-axis. Figure 3(b) shows the magnetic order is enhanced 
near the TBs. Defining the magnetic DWs to be where 
the magnetic order is suppressed, that is, near the middle 
between two TBs, then, clearly, the DWs are not located 
at the TBs. On the opposite sides of a DW, there is no 
change in the magnetic phase, thus we could label the 
DWs as the 0° DWs. Figure3(c) shows that the SC is en- 
hanced along the DWs, and that it is suppressed near the 
TBs. Figure 3(d) shows that the electron density is de- 
pleted near the TBs. Apparently, the depleted electron 
density leads to strong magnetic order that suppresses 
the SC order. On the DWs, the electron density appears 
to be close to optimal doping and thus SC gets enhanced. 

D. Diagonal twin boundary due to the asymmetry of 
As atoms- A TB due to missing one line of the lattice 
contains of both Fe and As (down) atoms oriented along 
45° from the x-axis is shown in Fig. 4(a). The D2d sym- 
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Figure 4: (a) The lattice structure near a twin boundary (red 
dashed line) with misplacing the As atoms along diagonal (or 
45°) direction. The open circles represent the positions of 
Fe atoms, and the red and blue dots respectively denote the 
As(up) and As(down) atoms. Corresponding spatial profiles 
of (b) the magnetic order, (c) the superconducting order, and 
(d) the charge density order are presented. 



metry of the lattice is also broken by the presence of this 
TB. Note that the geometry of this TB is fundamentally 
different from the one showed in Fig. 3 (a) since the TB 
does not pass through any of the Fe or As atoms. To 
study this case, we considered a 30x30 lattice with pe- 
riodic boundary conditions and three identical TBs ori- 
ented 45° from the x-axis (not shown). The TBs are 
located along y = x — 15, y = x and y = x + 15. 

Different from the case in Fig. 3(b), the magnetic or- 
der shown in Fig. 4(b) is suppressed along the TBs where 
DWs are located. The magnetic domain between the TBs 
still has the usual 2x1 collinear AF structure, except the 
magnetic moments are strongly and periodically modu- 
lated along the x-axis, which may be due to finite size 
effects of the TB and the distance between two near- 
est neighboring TBs. It also appears that the local 2x1 
collinear AF structure is replaced by a stripe-like y/2 x y/2 
AF structure oriented 45° from the x-axis. Furthermore, 
the SC is enhanced in these regions (see Fig.4(c)). From 
Fig. 4(d), note that on the TBs or DWs the carrier den- 
sity is corresponding to that in the overly hole-doped case 
(x ~ —0.3), which explains why the magnetic and SC 
orders are both suppressed on both sides of the TB. In- 
terestingly, stripe-like charge density waves oriented 45° 
from the x-axis occur on both sides of each TB. 

Our work has considered the effects of TBs on the com- 
plex interaction of magnetism and superconductivity in 
Fe-pnictides. There are three points that need to be em- 
phasized here. First, the formation and the location of 



the DWs are strongly related to the nature of TBs. For 
the four kinds of TBs studied here, the DWs in cases A, 
B and D are found to be pinned at the TBs, while in 
case C the DWs are separated from the TBs. Secondly, 
the formation of the DWs implies that the magnetism is 
inhomogeneous. This inhomogeneity strongly affects the 
SC and the electron density distribution. In cases of A, 
B and C, the SC is enhanced in the regions where m\ is 
suppressed and (m) is enhanced. The reasons for the en- 
hancement of SC are (i) the competition between the SC 
and magnetism, and (ii) the (n\) in the enhanced regions 
is close to the optimal doping level. However in case D, on 
the DWs where the m\ is suppressed, the electron density 
is close to over(hole)-doping level, which is unfavorable 
to SC. As a result, SC coexists with the magnetism in 
the middle of magnetic domains. Finally, we point out 
that our results on the diagonal TB of the lattice (se e 
Fig.l) are in good agreement with experiments (ToL [ill] . 
In a very recent scanning tunneling microscopy (STM) 
experiment for FeSe [20|, similar types of the TBs were 
detected but the SC was suppressed on the TBs. This 
may be caused by the differences in lattice structures 
and Fermi surface topologies of the two materials. We 
(l8| also find that the Coulomb-interaction-induced anti- 
phase DWs, as predicted in [8] and measured indir ectly 
in nuclear magnetic resonance (NMR) experiments [21] , 
always have higher ground energies than those of the 2x1 
collinear AF order and the 90° DW structures. These 
anti-phase DWs could form metastable states and may 
occur under certain conditions. In the present work we 
assume the lattices are well matched at the TBs. In the 
case that the lattice on both sides of the TB are not 
well matched, the SC could be suppressed by disorder 
scattering along the TBs. This issue deserves further in- 
vestigation. 
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